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ABSTRACT 



Context. High-resolution spectra of some chemically peculiar stars clearly demonstrate the presence of strong abundance gradients 
in their atmospheres. However, these inhomogeneities are usually ignored in the standard scheme of model atmosphere calculations, 
braking the consistency between model structure and spectroscopically derived abundance pattern. 

Aims. In this paper we present first empirical self-consistent stellar atmosphere model of roAp star HD 24712 with stratification 
of chemical elements included, and which is derived directly from the observed profiles of spectral lines without time-consuming 
simulations of physical mechanisms responsible for these anomalies. 

Methods. We used the LLmodels stellar model atmosphere code and DDAFIT minimization tool for analysis of chemical elements 
stratification and construction of self-consistent atmospheric model. Empirical determination of Pr and Nd stratification in the atmo- 
sphere of HD 24712 is based on NLTE line formation for Prii/iii and Ndii/iii with the use of the DETAIL code. 
Results. Based on iterative procedure of stratification analysis and subsequent re-calculation of model atmosphere structure we 
constructed a self-consistent model of HD 24712, i.e. the model which temperature-pressure structure is consistent with results of 
stratification analysis. It is shown that stratification of chemical elements leads to the considerable changes in model structure as 
to compare with non-stratified homogeneous case. We find that accumulation of REE elements allows for the inverse temperature 
gradient to be present in upper atmosphere of the star with the maximum temperature increase of about 600 K. 
Conclusions. 
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1. Introduction 

Atmospheres of chemically peculiar (CP) stars exhibit the 
presence of inhomogeneous distribution of chemical ele- 
ments. These inhomogeneities are frequently seen as abundance 
spots on stellar surfaces (see, for example, Kochukhov et al] 
l2004allbt iLueftinger et al.ll2003h and as clouds of chemical el- 
ements concentrated at certain depth s along the stellar radii 
dRvabchikova et al.ll2006i I2005L l2002h . Several tools have been 
developed and successfully applied to the analysis of element 
distributions based on detailed analysis of observed line pro- 
files. The horizontal elements distribution is investigated us- 
ing the rotational modulation of spectral lines: as star ro- 
tates, the abundance spots on its surface should reveal them- 
selves as variation of line profiles of respective elements ob- 
served at different rotational phases of the star The auto- 
matic procedure of finding the spots position and abundance 
values via the comparison of observed and calculated spec- 
tra is perf ormed via Dopple r Imaging techniques firstly intro- 
duced by [ Go ncharskv et al. (T982) and further developed by 
iPiskunov & Wehlau (1990) (see also Kochukhov et al.. (2004b)). 
The vertical distribution of chemical elements is modeled as- 
suming the stratification profiles of chemical elements to be 
represented by simple step-function with two different abun- 
dances in upper and lower atmosphere as was suggested in lBabell 
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(fT992h . or, as recendy showed bv lKochukhov et alj (l2006l) . ap- 
plying the regularized minimization algorithm. The results of 
these analyses show that strong abundance gradients do exist in 
atmospheres of CP stars with accumulation of various chemical 
elements at different atmospheric layers. 

From theoretical point of view, the appearance of inhomo- 
geneous element distribution with high abundance gradients in 
atmospheres of CP stars is explained by microscopic radiative 
diffusion processes: the atmospheres of these stars are believed 
to be dynamically stable enough for radiative levitation and grav- 
itational settling to rule the distribution of chemical elements 
(Michaud 1970). 

These abundance inhomogeneities may have rather strong 
influence on model atmosphere structure via modification of 
opacity and emissivity coefficients and thus the modeling of el- 
ement distributions in stellar atmospheres should, in principle, 
include the simulation of these diffusion processes on every step 
of model atmosphere calculation. Such models were success- 
fully developed and improved in last years (see, for example, 
Hui-Bon-Hoa et al. 2000; Alecian & Stift 2007), however, the 
sensitivity of modeled element distributions to the input physics 
(that is frequently poorly understood), the absence of accurate 
atomic data for some elements, and high computational expenses 
do not allow implementation of these models for detailed mod- 
eling of theoretical line profiles and their comparison with high- 
resolution observations: although the theoretical stratification is 
able to predict some observed characteristics of CP stars, it is 
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still not applied for the quantitative interpretations of modern 
high-resolution spectroscopic observations. 

One of the necessary ingredient in analysis of elements strat- 
ification is the computation of synthetic spectra and their subse- 
quent comparison with observations: a set of such computations 
is needed to construct the stratification profile of a given ele- 
ment that ensures a best fit between theoretical and observed 
line profiles. These empirical investigations of elements stratifi- 
cation in stellar atmospheres are very important for understand- 
ing the physical conditions in atmospheres of CP stars and, ide- 
ally, should provide the test ground for recent theoretical calcu- 
lations. The only methodological difficulty in empirical analysis 
is that almost in all cases the synthetic spectrum is computed 
using standard model atmospheres, i.e. model atmosphere that 
were computed under the assumption of homogeneous elements 
distribution. 

To avoid this inconsistency we made an attempt to imple- 
ment an iterative procedure of vertical stratification analysis with 
the subsequent re-calculation of model atmosphere structure. 
This procedure is applied to cool roAp star HD 247 12 for which 
the stratification of Si, Ca, Cr, Fe, Sr, Ba, Pr, and Nd is derived. 
Stratification analysis of Pr and Nd is based on non-local ther- 
modynamic equilibrium (NLTE) line formation for Prii/iii and 
Ndn/iii. Using this approach we examine how the stratification 
profiles of different elements are changing comparing with the 
case when only non-stratified model is applied and what is the 
overall effect of inhomogeneous elements distribution on model 
atmosphere structure. 

In the next section we give an overview of observations used. 
Then, in Sect. [3] we describe the general methods and atomic 
line data used for analysis of element stratification. In Sect. |4] 
we present the results of stratification analysis and its effect on 
model atmosphere structure and some observed parameters of 
the star The main conclusions are summarized in Sect. |5] 

2. Observations 

We used an average spectrum of HD 24712 obtained at 
November 10/11, 2004, with HARPS (High Accuracy Radial 
velocity Planet Searcher) spectrometer at the 3.6-m telescope 
at ESO, La Silla . The deta il s of t he reduction procedure are 
given in Rvabchi lcova et al.l (l2007h . In total 92 spectra were 
taken with 60s exposure time, signal to noise ratio S/N= 120, 
and R= 120000 resolving power These spectra cover a spec- 
tral region from 3850A to 6730A. After averaging the final 
S/N= 1000 was reached. HARPS spectra w ere obtained at the 
phase 0.867 close to magnetic maximum dRvabchikova et all 
l2007h . Photometric color-indices were taken from the SIMBAD 
elect ronic databas J| w ith the additional Stro mgren pho t ometr y 
from lMartine3 (Il993h . p ecuUar index a fromlVbgt et al.1 (1 19981) . 
and Geneva colors from B urki et al.l (11999 '). We also used UV 
spectroscopic observations extracted from lUE data archiv^. 

3. Methods 

3.1. Stratification analysis 

To analyse the stratification of chemical elements in the atmo- 
sphere of HD 24712 we applied step-function approximation as 
performed in DDAFIT - an automat ic procedure for de termina- 
tion of vertical abundance gradients (lKochukhovl2007h that was 

' http : // simbad . u- strasbg . f r/ simbad/ | 
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successfully used i n a number of studies (e.g. lRvabchikova et al.l 
120061 120051 120081) . In this routine, the vertical abundance dis- 
tribution of an element is described by four parameters: chemi- 
cal abundance in the upper atmosphere, abundance in deep lay- 
ers, the position of the abundance jump and the width of the 
transition region where chemical abundance changes between 
the two values. All four parameters can be optimized simulta- 
neously with the least-squares fitting procedure and based on 
observations of unlimited number of spectral regions. DDAFIT 
also allows the stratification analysis with the magnetically- 
splitted lines. For HD 24712 we assu med a pure radial magneti c 
field with the modulus <B)=3.1kG (Rvabchik ova et al.1 l2007h . 
The spectral synth esis was performed using SYNTHMAG code 
( lKochukhovl2007l) wh ich represents an improved version of the 
program developed bv iPiskunovl ( 1 1 9991) . 

For Pr and Nd, we performed NLTE str atification analysis 
as described in iMashonkina et al] ( l2009l) and lMashonkina et al.l 
12005) using a trial-and-error method and the observed equiva- 
lent widths of the lines of the first and second ions. 



3.2. NLTE calculations 

The present investigation of Pr n/m and Nd n /m is based on the 
NLTE met hods treated in our earlie r studies dMashonkina et al.) 
200^ and (IMashonkina et alj|2005h . where atomic data and the 
problems of line formation have been considered in detail. To 
describe briefly, comprehensive model atoms of praseodymium 
and neodymium include the measured and the predicted energy 
levels, in total, 6708 levels of Prii and Priii and 2258 levels 
of Nd II and Nd iii. The coupled radiative transfer and statisti- 
cal equilibrium equation s were solved using a revised version 
of the DETAIL program dButler & Giddingsll l98 5) based on the 
accelerated l ambda iteration, which follows the efficient method 
described by'Rybicki & Hummeil ([19911 ll992|). 

As shown for Nd by Masho nkina et alj (l2005h and for Pr by 
IMashonkina etalJ(l2009h .the main non-LTE effect for their first 
ions in the model atmosphere representing the atmosphere of 
HD 247 12, is overionization caused by a super-thermal radiation 
of non-local origin close to the thresholds of the low-excitation 
levels. The departures from LTE for the lines of the first and 
the second ions are of the opposite sign and they are signifi- 
cant if the element is concentrated in the uppermost atmospheric 
layers where collisions are inefficient in establishing thermody- 
namic equilibrium. For the Pr and N d distribution in the atmo- 
sph ere of HD 24712 determined by IMashonkina et all d2009l) 
and lMashonkina et alj d2005l) . the NLTE abundance corrections, 
AeNLTE = ^NLTE - ^LTE, are positive for the lines of Prii and 
Nd II and constitute between 1 .0 dex and 1 .4 dex, while AeNLXE 
is negative for the lines of Pr iii and Nd m and ranges between 
-0.3 dex and -0.7 dex. 

We perf ormed NLTE calculation s for Hi using the method 
described in lMashonkina et all (12008). The model atom includes 
levels with principal quantum numbers up to n - 19. For 
the model atmospheres investigated in this study, the ground 
state keeps their thermodynamic equilibrium level populations 
throughout the atmosphere. In the layers where the core-to-wing 
transition is formed, namely, between logrgooo - -1-2 and 
log T5000 - -2.2, the n - 2 level is closely coupled to the ground 
state, while departures from LTE for the n - 3 level are con- 
trolled by Ha which serves as the pumping transition resulting 
in an overpopulation of the upper level. For Ha, NLTE leads to 
weakening of the core-to-wing transition compared to the LTE 
case. 
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3.3. Model atmospheres 

To perform the model atmos phere calculations w e used the re- 
cent version of the LLmodels dShulvak et al.l200 4^ stellar model 
atmosphere code (version 8.6). LLmodels is 1-D, hydrostatic, 
plain-parallel LTE code which accounts for the effects of individ- 
ual and stratified abundances. It is to note that the stratification 
of chemical elements is an input parameter for LLmodels code 
and thus is not changing during model atmosphere calculation 
process. This allows to explore the changes in model structure 
due to stratification that was extracted directly from observations 
without modeling the processes that could be responsible for the 
observed inhomogeneities. 

The following calculation settings were used for every model 
atmosphere calculation: the atmosphere of a star is discretized 
into 120 layers between logrgooo - -8 and 2 with a non- 
equidistant spacing, i.e. with the higher points density in the re- 
gions of steep abundance gradients to ensure accurate integra- 
tion of radiation field properties and solution of other model 
equations. The frequency-dependent quantities are calculated 
with variable wavelength step: 0.1 between 500A and 20000A, 
and 0.5A between 20000A a nd 70000A respec tively. VALD 
database (Piskunov et al. 1995t iKupka et al.|[l999l) was used as 
a source of atomic hnes data for computation of Une absorption 
coefficient. 

The use of the LTE assumption in atmosphere modeling re- 
sults from the following consideration. The opacity in the at- 
mosphere of HD 24712 is mostly contributed by neutral hydro- 
gen. Our NLTE calculations for Hi showed no significant de- 
viations from thermodynamic equilibrium populations for the 
n — \ level throughout the atmosphere and for the n - 2 
level below log 75000 = -3. The departures from LTE are found 
to be larger for the hotter at mosphere of Vega , however, self- 
consistent NLTE modelling of iHauschildt et all (Il999h leads to 
only 200 - 300 K changes in T distribution compared to the LTE 
model in the layers between logrsooo = ^1 and logrsooo - -6. 
The NLTE effect on atmospheric structure of HD 24712 is ex- 
pected to be smaller. Another important opacity contributors are 
the iron group elements. The iron and chromium among them 
are concentrated in deep atmospheric layers of HD 24712, be- 
low logTsooo - -2 (see Fig. [3]), where the departures from LTE 
are expected to be small. This is confirmed by our NLTE cal- 
culations for Ca which distribution in the atmosphere is similar 
to that for the iron group eleme nts (Fig.[3l). We used th e NLTE 
method for Ca i/ii described by Mashontdna et alj (|2007|) . We do 
not know the dist ribution of light elemen ts, C, N, and O, how- 
ever, according to lKhan & Shulvafl (l2006h these elements do not 
introduce noticeable changes in model structure and emergent 
flux. The influence of the departures from LTE for two rare- 
earth elements (REE), Pr and Nd, on atmospheric structure of 
HD 24712 is simulated in present paper. 

Finally, the magnetic field with the intensity less than 5 kG 
has marginal impact on model temperature-pressure structure 
d^chukhov et al. 2005; Khan & Shulyak 2006). Thus, we ig- 
nored the effects of polarized radiative transfer and Zeeman 
splitting in all model atmosphere calculations presented in this 
study. 

3.4. Atomic line data for stratification analysis 

As was emphasized in 'Rvabchikov a et al] (l2003h the vertical 
abundance stratification manifests itself as an impossibility to fit 
the core and wings of strong lines with developed Stark wings 
(Can K, Sin, Mgii lines, etc.) with the same abundance or as 
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Fig. 1. Fe abundance distribution derived with the homogeneous 
model atmosphere and the range of the optical depth formation 
of Fe I and Fe 11 lines from Table[T]calculated with homogeneous 
Fe abundance log(Fe/A^totai) - -5.10. 

an impossibility to describe the low and high-excitation lines 
of the same ion with a chemically homogeneous atmosphere. 
It means that the choice of a proper set of spectral lines is cru- 
cial for element stratification study. Ideally, it should be a set 
of lines originated from the levels in large range of the excita- 
tion energies and from different ions, thus probing a substantial 
part of stellar atmosphere. In chemically peculiar stars as cool 
as HD 24712 we deal with severe blending by numerous lines 
of the rare-earth elements, therefore the strongest lines with the 
developed wings (like, for example, Fe i A 4045.8 A line) are not 
possible to use. Also, Fe 11 lines with the upper level excitation 
potential at 10 eV and higher are not visible any longer The 
blue spectral region is more crowded with the REE lines than 
the red one and the preferential working spectral region is 5000- 
6800AA. Table[T]gives a list of Si, Ca, Cr, Fe, Sr, Ba Unes used in 
LTE stratification analysis routines together with the atomic pa- 
rameters and the references. For NLTE stratific ation analysis we 
used t he sa me lines as given i n the p apers of .Mashonkina et all 
(120051) and lMashonkina et al.l (l2009b . 

An ability of the chosen set of spectral lines to probe dif- 
ferent atmospheric layers is illustrated in Fig.[T] The iron strat- 
ification derived with the initial model (see Section|4]i is shown 
together with the ranges of Fe lines depth formation calculated 
on log Tgooo scale. The later were calculated with t he contribution 
function to the emergent line radiation following 'Achma d et al.l 
(1991). Fig. |2] illustrates a fit of the calculated to the observed 
Fe Line profiles for homogeneous and stratified atmospheres. For 
Si, Ca, Cr, Sr, and Ba the corresponding plots are shown in 
Figs. |7|8l9|10l (Online material). 

3.5. Self-consistent models with empirical stratification 

Using the methods of stratification analysis described above one 
could restore the element distribution profile of any chemical 
elements for which accurate atomic line data exist. However, 
we should remember that the empirical analysis of chemical el- 
ements stratification is based on model atmosphere technique, 
and thus the temperature-pressure structure of model atmo- 
sphere itself depends upon stratification which has to be found. 
Therefore, the calculation of model atmosphere and stratification 
(abundances) analysis are linked together and thus an iterative 
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Table 1. A list of spectral lines used for the stratification calculations. The columns give the ion identification, central wavelength, 
the excitation potential (in eV) of the lower level, oscillator strength (log(g/)), the Stark damping constant ("appr" marks lines for 
which the classical approximation was used), and the reference for oscillator strength. 
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G - 'Garz (1973); BBC - Blanco et al. (1995); SG - Smith & Gallagher (1966); SMP - Seaton et al.' ('1994')); S - Smith (1988)); SO - 
[Smith & O'Neil ( 1975); SR - Smith & Raggett ( 1981); RU - Raassen & Uvlings ( 1998); MWF - Martin et al. ( 1988); NL - Nilsson et al ( 20(}d ); 
BWL-|Black et al. ( 1972); BPS 1 - Blackwell et al. ( 1979 ); BPS2 -LB lackwell et al. ( 1982); BK-Bard & Kock ( 199 4i); BKK-lBard et alJ J1991) : 
K07 - (Kurucz databastO) W - .Warner (1968) ; GC - .Garcia & Campos (1988) ; * - .Fleurier et al., ( .1977,) ; NBS - .Miles & Wiesd ( Il969l) . 



procedure should be used in this case. In general, this procedure 
could be divided into several steps: 

1. calculation of so-called first approximation, simplified 
model, that could be any standard model atmosphere with 
the atmospheric parameters close to those of the star to be 
analysed; 

2. determination of stratification of chemical elements based on 
spectroscopic analysis of individual spectral lines; 

3. calculation of improved model atmosphere taking into ac- 
count stratification found in previous step; 

4. comparison of modeled energy distribution (or/and photo- 
metric colors) and hydrogen line profiles with observed ones 
and re-adjustments of model input parameters {T^sf, \ogg) if 
needed to fit observations. At this step several calculations 
of model atmosphere with fixed stratification but different 
model parameters may be required; 

5. repeating the overall process starting from step 2 until strati- 
fication profiles of chemical elements and model parameters 
are not converged. 



At the end, the listed above procedure ensures the consis- 
tency between model atmosphere structure and abundances used 
for calculation of synthetic line profiles and interpretation of ob- 
served data. 

4. Results 

4.1. Construction of model atmosphere 

In the present analysis we started from chemically homoge- 
neous model atmosphere with the parameter s Teff - 7250 K and 
log^ - 4.3 and the abundances taken from 'Ryabc hikova et alJ 
(1997) (Model 1). These mean abundances were derived based 
on the observations close to magnetic maximum, and they pro- 
vide a satisfactory fit for HARPS observations. Table |2] shows 
abundances of the elements which were not passed through the 
subsequent stratification analysis. For the rest of chemical ele- 
ments the solar composition ( Asplund et al. 2005) was assumed. 

Based on this model the stratification of Fe, Ba, Ca, Sr, Cr, 
Si, Pr, and Nd was derived and is displayed on the top panel of 
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Wavelength 

Fig. 2. A comparison between the observed Fe line profiles and calculations with the stratified abundance distribution (full line) and 
with the homogeneous (log(Fe/A^totai) = -5.10) abundances (dashed line). 

Fig. |3] It is seen that the stratification of majority of elements 
shows element underabundance in surface layers and overabun- 
dance in layers around star's photosphere (logTsooo ~ [-1,0]). 
Contrary, the distribution of Pr and Nd show inverse picture 
with the strong overabundance (up to 5dex) in surface lay- 
ers. This kind of elements distribution, i.e. the difference be- 
tween distribution of REE's comparing to other elements, seems 
to b e a common characteristic of the atmospheres of Ap stars 
jRyabchikova et al. 2002). 

The stratification shown in Fig. |3] (top panel) was then ap- 
plied to the calculation of new model atmosphere. Again, the 
individual values of abundances of some other elements used for 
calculation were taken from Table |2] 

The stratification of REE elements and its implementation in 
model atmosphere computation deserve more detailed descrip- 
tion. Numerical calculations showed that the REE elements strat- 
ification has a strong influence of model atmosphere structure 
leading to the appearance of inverse temperature gradient in up- 
per atmospheric region. This is illustrated in Fig. |4] where the 
models with the same atmospheric parameters and input line list 
were computed with and without stratification of chemical el- 
ements (solid and dotted lines respectively). It is seen that the 
heating of the stratified atmosphere starts at log T5000 - -3 going 
outwards to the stellar surface and reaches its maximal of about 
800 K around log tsooo = -4 compared to homogeneous abun- 
dance model. This temperature increase is caused by the pres- 
ence of REE stratification and, as it seen from Fig. [3] is directly 
related to the position of abundance jumps of Pr and Nd. No 
other elements considered in our computations are responsible 
for such an effect. Indeed, the strong overabundance of Pr and 
Nd by 4-5 dex relative to their solar values and their rich spectra 
presented in current version of VALD database make the absorp- 



Table 2. Abundances of individual elements. 



Element log(jVei/Motai) log(A^ei/Motai)o 



c 


-4.00 


-3.65 


N 


-4.50 


-4.26 





-4.00 


-3.38 


Mg 


-5.60 


-4.51 


Al 


-5.57 


-5.67 


Ti 


-7.28 


-7.14 


Mn 


-7.01 


-6.65 


Co 


-5.57 


-7.12 


Ni 


-6.54 


-5.81 


Y 


-7.80 


-9.83 


La 


-9.00 


-10.91 


Ce 


-8.90 


-10.46 


Sm 


-9.00 


-11.03 


Eu 


-9.30 


-11.52 


Gd 


-8.70 


-10.92 


Dy 


-8.94 


-10.90 


Er 


-9.53 


-11.11 



The ab undances derived from observations are from lRvabchikova et al.l 
1 I1997I) . The REE's abundances were deriv ed using the lines of the first 
ions. The solar abundances are taken from lAsplund et al.1 d2005h . 



tion coefficient to be very high at frequencies of REE transitions. 
The latter are generally located at the spectral region where the 
star with Tgff - 7250 K radiates most of its energy dominating 
in the total radiative energy balance, and thus allowing for the 
heating of atmospheric layers with high REE opacity. 

The appearance of inverse temperature gradient clearly 
shows that the REE elements are no longer trace elements in at- 
mospheres of Ap stars if their stratification is as strong (or com- 
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Fig. 3. Stratification of eight elements in atmosphere of 
HD 24712 derived using homogeneous first-approximation 
model (Model 1, top panel) and using final model with scaled 
REE opacity (Model 5, bottom panel). 
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Fig. 4. Effect of Pr and Nd stratification on model temperature 
structure. Solid line - t7250g4.3 model calculated with homo- 
geneous abundances; dotted line - t7250g4.3 model calculated 
with stratified abundances shown in Fig. |3j dashed and dash- 
dotted lines - t7250g4. 1 models calculated with stratified abun- 
dances and scaled Pr n,m and Nd ii,m logigf) values. 



parably strong) as in the atmosphere of HD 24712. This immedi- 
ately implies that not only synthetic line profiles, but also model 
atmospheres with REE stratification must be computed based on 
NLTE line formation for REE's. Unfortunately, this can not be 
done with the current version of LLmodels. Nevertheless, to sim- 
ulate NLTE effects of Pr and Nd in model atmosphere calcula- 
tions we used the following approach. iMashonkina et all (l2005h 
and Mashonkina et al. (2009) found that, for the stratified distri- 
bution of Nd and Pr, the NLTE abundance corrections are posi- 
tive for the lines of the first ions at a level of one order of mag- 
nitude and negative for the lines of the second ions at a level of 
^ 0.5 dex. Thus, to account for the REE stratification in LTE 
model atmosphere code, we scale the log(,g/) values of Pr ii and 
Ndii lines by -L5dex while taking abundances derived using 
second ions for model atmosphere calculation. This procedure 
was applied to all the Pr ii and Nd ii lines that are presented in 
master line list used for model atmosphere calculations. Then, 
using this new line list, we recomputed the model atmosphere 
and re-derived stratification of chemical elements again. At each 
iteration Pr and Nd stratifications were recalculated based on 
NLTE line formation for Pr ii/iii and Nd ii/iii. The best fit to hy- 
drogen Ha line profile (wings) and photometric color-indices 
required the decrease of effective gravity down to \ogg = 4.1 
(mainly to fit c\ index). No further noticeable changes were de- 
tected in observed parameters and this gravity was considered 
as a final one. The model and stratification calculations are con- 
verged after the second iteration (Model 4). 

Except Pr and Nd, other REE elements can also be non- 
uniformly distributed in the atmosphere of HD 24712 and can 
also exhibit deviations from LTE. Among them, only most abun- 
dant elements may play a noticeable role in modification of at- 
mospheric structure. At present, we do know neither distribu- 
tions of these elements, nor the magnitude of their NLTE effects. 
In any case, preliminary abundance results for Ce and Sm show 
that REE anomaly - a discrepancy between abundances inferred 
from the lines of singly- and doubly -ionized element - reaches 
approximately the same magnitude as for Pr and Nd. Supposing 
that the Ce (which has the biggest number of available lines in 
VALD among other REE elements) has the same vertical distri- 
bution as Nd, we estimated its effect on atmospheric structure. 
We found that neglecting Ce stratification we underestimate a 
size of temperature jump by no more than ~200 K. 

Decreasing the log(g/) values of Prii and Ndii lines could, 
in principle, simulate the NLTE line opacity of these elements, 
however the second ions are also affected by NLTE effects that 
bring as much as 0.5 dex difference between LTE and NLTE 
abundances derived using second ions spectra. This forced us 
to compute another model where the log(^/) values of both 
first and second ions are changed by -Idex and 0.5 dex re- 
spectively (Model 5). This model provides slightly higher tem- 
perature increase in superficial layers compared with the model 
where log(,g/) values of only first ions were changed (see Fig.|4]l. 
It is seen that the simulation of NLTE REE opacity tends to de- 
crease the amplitude of the temperature jump almost by a factor 
of two. However, independently on the way how to scale the 
REE opacity, the overall picture remains the same illustrating 
the heating of surface layers in the region of strong Pr and Nd 
overabundance. 

We also tested the influence of Prii and Ndii bound-free 
opacity on model structure and energy distribution. To simplify 
the computations we assumed three characteristic energy levels 
with excitation energies of 2 eV, 3 eV, and 4 eV for both ions. 
These levels further consist of 50 atomic states each with the sta- 
tistical weight g -1 ■ Such an approximation is based on the fact 
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Table 3. Observed and calculated photometric parameters of HD 24712. 







t7250g4.3 


t7250g4.3 


t7250g4.1 


t7250g4.1 


t7250g4.1 


Color 




homogeneous 


stratified 


stratified 


stratified 


stratified 


index 


Observations 




REE unsealed 


REE unsealed 


Alog(g/)p,.n,Ndii = -1.5 


Alog(g/)PiII,NdII = -1 

Alog(g/)p,.m.Ndiii = +0.5 






(Model 1) 


(Model 2) 


(Model 3) 


(Model 4) 


(Model 5) 


b-y 


0.186 (0.003); 0.191* 


0.195 


0.188 


0.183 


0.180 


0.180 


m\ 


0.202(0.006); 0.211* 


0.197 


0.214 


0.213 


0.205 


0.207 


Cl 


0.653(0.011); 0.626* 


0.610 


0.575 


0.638 


0.663 


0.656 


Hp 


2.745(0.003); 2.760* 


2.802 


2.804 


2.808 


2.810 


2.809 


a 


0.609 (0.0005) 


0.621 


0.621 


0.621 


0.617 


0.619 


B-V 


0.320 


0.308 


0.308 


0.297 


0.288 


0.290 


U-B 


1.381 


1.333 


1.311 


1.370 


1.379 


1.376 


V-B 


0.572 


0.596 


0.595 


0.610 


0.621 


0.619 


Bi-B 


0.978 


0.984 


0.990 


0.985 


0.982 


0.983 


B2-B 


1.393 


1.413 


1.407 


1.413 


1.417 


1.416 


VrB 


1.286 


1.319 


1.319 


1.333 


1.343 


1.341 


G-B 


1.688 


1.727 


1.726 


1.741 


1.753 


1.751 



* - lMartinezlfl993h 

The values in brackets give the error bars of observations. 



that both ions have essentially the same structure of states with 
the very close ionization potentials. The NLTE level populations 
were taken directly from DETAIL code and then incorporated in 
LLmodels. For all the levels, the photoionization cross-sections 
were calculated using the hydrogenic approximation. It appeared 
that the effect of the Pr 11 and Nd 11 bound-free opacity becomes 
visible only if the abundances of these elements are increased 
by -1-6 dex for the stratified model where REE's are already in 
strong overabundance in surface layers. Even so, the cumulative 
effect on on the temperature distribution was never higher than 
tens of K. We, thus, concluded that REE bound-free opacity can 
be ignored in the model atmosphere computations. 

Bottom panel of Fig. [3] shows the stratification derived with 
Model 5. One can note that the stratification of REE's shows 
more steep abundance jump in this case which is probably a re- 
sult of more refined set of optical depths around REE jump im- 
plemented in the last iteration compared to the initial model and 
the first iteration. This is also the reason why the inverse temper- 
ature gradient shown in Fig.|4]is more steep for the model calcu- 
lated with Model 5 stratification than for the model with strati- 
fication derived in Model 4. Few remarks concerning Sr and Ba 
stratification calculations should be made. Although we get for- 
mally Ba stratification, it is not significant because the difference 
between for line profiles fits with uniform and stratified mod- 
els does not exceed 15-16%, while differs by 5-10 times in 
the case of real stratification. Therefore uniform Ba distribution 
in the atmosphere of HD 247 12 or at least in the layers where the 
studied Ban lines are formed is more probable. For Sr distribu- 
tion we used five lines with three weakest lines of Sri. Strontium 
distribution is mainly based on strong Sr 11 lines, in particular, on 
resonance Srii A 4215.5 A. The wings of Srii lines that probe 
abundance in the deeper atmospheric layers, are blended, that 
obviously leads to larger formal computing errors of stratifica- 
tion parameters determination and, hence, to worse To check 
a reliability of Sr stratification we recalculated it using Sr lines 
observed at slightly differe nt magnetic phase . 944 (s ee UVES 
time-series observations in iRvabchikova et al] (l2007h ). Within 
the formal errors of DDAFIT procedure strontium abundance 
profiles coincide below log T5000 = -2 where practically all con- 
sidered Sr lines are formed. They differ above log rsooo - -2, 
where the only core of Srii A 4215.5 A line is formed. Poor fit of 



this line in the transition region between the line core and the line 
wings indicates that a more complex shape of abundance profile 
than a simple step-function is required to desc ribe the observed 
line profile, similar to the Ap star HD 133792 (iKochukhov et al.l 
2006). Strontium stratification derived for HD 24712 is not spu- 
rious, but the shape of the Sr abundance profile may differ from 
that derived in the present study. The other elements that demon- 
strated strong changes during iterative stratification analysis are 
Si and Ca. Both elements show the increase of the abundance 
jump amplitude in the final stratification with strong underabun- 
dance in surface layers and with practically the same abundance 
near photospheric layers. The jump itself becomes steeper. It 
may be caused by the change in temperature and electron den- 
sity structure of the atmosphere above log tsooo = -3, where the 
central parts of the strongest Cai A 4227 and Sin AA 6347,6371 
are formed. However, the formal errors of Si and Ca abundances 
in the upper atmospheric layers derived in DDAFIT procedure 
are rather high, up to 1.5 dex, hence 2.5 dex difference for the 
initial and final Ca distributions may be considered as within the 
3cr limit. 

The changes in element distributions between the two cases 
shown in Figs. [3] are purely due to the changes in model atmo- 
sphere structure. Thus we can conclude that the recalculation of 
the model atmosphere structure could be very important for the 
quantitative stratification analysis. For some elements, the strat- 
ification derived with models computed under the assumption of 
homogeneous elements distribution can be considered only as a 
first approximation. 

Finally, in Table[3]we compare theoretical and observed pho- 
tometric color-indices of HD 24712. The calculations are pre- 
sented for two models with T^ff = 7250 K, logg - 4.3 (homo- 
geneous and stratified abundances) and for three models with 
Teff - 7250 K, \ogg = 4.1 but with different assumptions 
about Pr and Nd opacity. The theoretical co lors were ca lculated 
with modified computer codes taken from Kuruczi (Il993a ). One 
can see that the color-indices of Stromgren system as well as 
U-B and V-B indices of Geneva system are better repre- 
sented by stratified models . In contrast, the last two Geneva 
indices given in Table [3] are better described by homogeneous 
abundance model. However, these differences are not very big 
and we can conclude that, generally, the two final models cal- 
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culated with the complex stratification shown in Fig. |3] (bottom 
panel) and scaled REE opacity demonstrates a good agreement 
with the observatio ns. We have to no te that the Stromgren pho- 
tometry taken from iMartined (1 19931) . which is not available in 
electronic format, int roduces addition al uncertainties. In partic- 
ular, index b - y of lMartine3 (Il993h is better represented by 
homogeneous model. Furthermore, the difference of 0.027 mag 
between the two observations for ci index (see Table O does 
not allow to distinguish definitely between homogeneous and 
strati fied models. On the other ha nd, the observation of ci taken 
from lHauck & MermilliodI (Il998l) agrees nicely with predictions 
by model with scaled Prii,iii and Ndii,iii opacity. The differ- 
ence of 0.007 mag between models computed under two as- 
sumptions about the scaling of REE opacity (last two models 
in Table O is within the error bars of observations, but devi- 
ates significantly from ci calculated without REE scaling. Our 
calculations demonstrate that taken stratification effects into ac- 
count we obtain lower surface gravity than follows from the 
cahbrations based on homogeneous models. Strong impact of 
simulated NLTE effects of Pr and Nd on ci value, tempera- 
ture distribution, and colors demonstrates a necessity of more 
precise model atmosphere computations with accurate treatment 
of Pr, Nd, and possibly other REE lines formation. Our choice 
of logg for the final model could also change as more ad- 
vanced NLTE model atmospheres are employed. Moreover, the 
strong variabi lity in ci index with the peak-to-peak amplitude 
of 0.086 mag dWoUf & MorrisonllT973l) is likely to be resp onsi- 
ble for the difference betwe en the value of iMartinez* ('1993|) and 
the average value given bv iHauck & Merm illiod (1998). Thus, 
at present it is difficult to say something more definite about the 
log g value of the star and we restrict ourselves with log g =4.1. 

As it was stated above, for the opacity calculation in model 
atmospheres we used REE data extracted from VALD database. 
However, for the NLTE calculations with DETAIL code it 
was necessary to have more extensive line lists for accurate 
calculations of radiative rates for every single atomic state. 
These detailed line lists for Prii and Ndii/iii were produced 
in Institute of Spe ctroscopy (ISAN, Russia), and described in 
iMashonkina et al.1 (2005, 2009). Note that the ISAN data al- 
ready contains Pr and Nd parameters from VALD thus providing 
only additional number of predicted lines. For example, while 
VALD provides information about 508 Prii, 1279 Ndii, and 55 
Ndm lines, ISAN theoretical computations extended them up 
to 103428, 1172579, and 6858 lines respectively. However, this 
great increase in line numbers showed only marginal changes 
in model temperature distribution. This allowed to conclude that 
the radiative equilibrium balance in atmosphere of HD 24712 is 
controlled by strongest REE lines that are presented in VALD. 
The only noticeable changes were detected for the ci (increased 
by 0.012 mag) and for mi index (decreased by 0.008 mag). 

Fig. |6] demonstrates the observed and predicted energy dis- 
tributions of HD 24712. Observational data is represented by 
a combination of the lUE low-resolution flux and photometric 
measurements converted to absolute flux units. The theoretical 
fluxes are calculated with the last model from Table[3] Note that 
all the models from Table[3]represent observations good enough, 
and the only possibility to choose the best one would be the com- 
parison of U V fluxes blueward A = 2000A where stratified mod- 
els predict significant flux excess (almost one order of magnitude 
in maximum) compared to homogeneous abundance model, but 
this region plays unimportant role in the total energy balance in 
the atmosphere due to the relatively cool effective temperature of 
the star Although the lUE observations tend to be closer to the 
theoretical stratified fluxes at this region, it is difficult to make 
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Fig. 5. Observed and calculated Ha line profiles. Thick full 
line - HARPS observations; dotted line - model with Teff = 
7250 K, log,g = 4.3 and with homogeneous elements distribu- 
tion; full line - model with Teff = 7250 K, log^ = 4.3 and with 
stratification shown in|3](top panel); dashed line - model with 
Teff = 7250 K, log^ = 4.1 and with final iteration shown in|3] 
(bottom panel). Bottom panel shows the zoomed part of the fig- 
ure on top panel around the core of the line. 



a clear decision because of large error bars in observations (see 
lUE data archive). 

Fitting of the spectral energy distribution of HD 24712 cal- 
ibrated in physical flux units allows us to derive the angular di- 
ameter of the star, = 0.335+ 0.005 mas. Combin ed with the re- 
vised Hipparcos parallax of tt = 20.32 + 0.39 mas dvan LeeuwenI 
'2007|), this yields a new determination of the stellar radius, 
1.772 + 0.043/?o- 

4.2. Effect of REE stratification 

It was shown above that the accumulation of Pr and Nd leads to 
the appearance of inverse temperature gradient with the temper- 
ature jump of about 600 K as compared to non-stratified model. 
However, this temperature jump is located in upper atmospheric 
layers and has no influence on the wings of Ha line, as it is 
shown in Fig.|5] Thus, even in the case of strong stratification of 
REE elements, the determination of atmospheric stellar param- 
eters by fitting the wings of hydrogen lines could be performed 
with the model not taking the stratification of REE into account. 
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Fig. 6. Calculated and observed energy distributions of HD 24712. Solid line - t7250g4.1 model calculated with stratified abun- 
dances shown in Fig. [3] (bottom panel) and changed Pr ii,iii and Nd ii,iii log(g f) values by - 1 dex and +0.5 dex respectively. 



However, LTE calculations show that the core of Ho' line is 
influenced strongly by stratification of Pr and Nd (bottom panel 
of Fig.|5]l. One could see that none of the He profiles computed 
under the assumption of LTE line formation could fit the transi- 
tion region between the line core and wing, that is commonly 
called as core-to-wing anoma ly (CWA) and is freque ntly ob- 
served in spectra of Ap stars. iKochukhov et alj (|2002|) tried to 
explain this anomaly by the existence of the temperature bump 
in the upper atmospheric layers above logrsooo ~ -1- What is 
interesting is that the heating of atmospheric layers due to REE 
stratification produces the same sharp transition between the line 
wings and core, but the central intensity of theoretical spectra are 
much shallower than that with homogeneous abundance model 
which still provides the best fit. However, all stratification ef- 
fects disappear when NLTE hydrogen lines formation is applied, 
which is 'washing out' all stratification effects. 

From Figs. 2 and 3 of IKochukhov et aTl (|2002|) it is seen 
that in order to explain the core-to-wing anomaly of hydrogen 
Balmer lines, the temperature jump should be located between 
logTsooo ~ -3 and -2. In our calculations the temperature gradi- 
ent starts higher in stellar atmosphere at log tsqoo k -3. Position 
of the temperature gradient caused by REE stratification depends 
on the details of NLTE line formation, but it is practically im- 
possible to shift REE jump position in HD 24712 downward 
by m ore than 0.3 - 0.5 dex in log T5000 scale dMashonkina et alj 
120091) . 

REE stratification is not the only reason for inverse temper- 
ature g radie nts in Ap star atmosphere. Recently, Alecian & Stift 
(I2OO7I) and iLeBlanc et alj (l2009i) calculated element stratifica- 
tion in the atmospheres with magnetic field and showed that hor- 
izontal magnetic field changes abundance distribution. Instead 
of a step-like shape which is obtained in diffusion calculations 
in non-magnetic atmospheres or in the presence of predomi- 
nantly vertical component of the field, one gets an abundance 
minimum and then an increase of element abundance in the up- 



permost atmospheric layers. These abundance increase might 
lead to the corresponding te mperature increase sta rted just at 
logTsooo = -1 (see Fig. 5 of iLeBlanc et al.l l2009l) . In a sim- 
ple step function approximation used in the present s tratification 
analysis as well as in the most other similar studies dWade et all 
2001; Ryabchikova et al. 2005, 2006, 2008), we cannot probe 
an increase of element abundance in the uppermost atmospheric 
layers. However, relatively good fitting of the observed line pro- 
files in many stars with different strengths of magnetic field ob- 
tained with step-like abundance distribution of Si, Ca, Cr, Fe 
- main contributors to line opacities - supports the absence of 
significant increase of element abundance immediately follow- 
ing the abundance jump. This is further corroborated by the 
vertical inverse problem solution for the Ap star HD 133792 
(IKochukhov et alj |2006), where the authors were able to derive 
element distribution without a priori assumption about the shape 
of this distribution. One should note, however, that this method 
can be currently applied only to non-magnetic or weakly mag- 
netic atmospheres. Additionally, in atmospheric diffusion calcu- 
lations the stratification is generally obtained assuming that all 
the forces acting on a given plasma volume are in stationary 
equilibrium (total diffusion velocity is zero). Regarding to real 
stars this is not necessarily the case and element stratifications at 
a given moment could be different (see lAlecian & Stift 2007). 

Thus, to explain the most part of the observed peculiarities 
in cool Ap stars we need to implement the full NLTE treatment 
of REE elements in model atmosphere calculations and not the 
rough simulation presented in this study. The accumulation of 
other elements could also be responsible for CWA and more de- 
tailed investigation is needed to explore this effect. 

5. Conclusions 

The spectroscopic investigation of elements stratification is the 
only method for testing modern predictions of the particle dif- 
fusion theory in atmospheres of star where this mechanism is 
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sufficient to produce detectable abundance gradients. It provides 
us with stratification profiles of different chemical elements that 
could be compared with recent self-consistent diffusion models. 
However, in most of the stratification analysis routines the stan- 
dard stellar model atmospheres are used where the homogeneous 
abundances are assumed. 

To circumvent this difficulty we made an attempt to construct 
an empirical self-consistent model atmosphere where the strati- 
fication of chemical elements is derived directly from observed 
spectra and then treated in a model atmosphere code. This itera- 
tive procedure was applied to one of the cool CP star HD 24712. 
Below we summarise the main results of this investigation: 

- The accumulation of Pr and Nd in the upper atmosphere has 
an impact on the temperature structure of the atmosphere. 
The strong overabundance of these elements leads to an ap- 
pearance of inverse temperature gradient with the maximum 
temperature increase of up to 600 - 800 K as compared to 
homogeneous abundance model. 

- For HD 24712 we find that the effect of elements stratifica- 
tion has minor influence on the observed energy distribution. 
The lack of good quality observations in UV region (where 
the effect of stratification is more pronounced) does not al- 
low at present to carry out more precise investigations. The 
changes in temperature structure due to the stratification of 
Pr and Nd do not affect the wings of hydrogen lines and thus 
can be ignored in routine procedures of the determination of 
atmospheric stellar parameters. 

- Due to their high overabundance and important role in over- 
all temperature balance, REE element can no longer be con- 
sidered as trace elements in standard scheme of NLTE cal- 
culations. This implies that, to achieve a better consistency 
in modelling of REE stratification, the model atmosphere 
should account for REE NLTE effects as well. 

- For HD 24712 we find no critical changes in the shape 
and amplitude of stratification profiles of chemical elements 
when recalculating the atmospheric temperature-pressure 
structure. Most of the changes are comparable with the er- 
ror bars of the minimization procedure used. 
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Fig. 7. A comparison between the observed Si Line profiles and calculations with the stratified abundance distribution (full line) and 
with the homogeneous (log(5 //A^totai) - -4.50) abundances (dashed line). 
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Fig. 8. A comparison between the observed Ca Line profiles and calculations with the stratified abundance distribution (full line) and 
with the homogeneous (log(Cfl/Motai) = -5.80) abundances (dashed line). 
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Fig. 9. A comparison between the observed Cr Line profiles and calculations with the stratified abundance distribution (full line) and 
with the homogeneous (log(Cr/A^totai) - -5.80) abundances (dashed hne). 
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Fig. 10. A comparison between the observed Sr(left panel) and Ba (right panel) line profiles and calculations with the stratified 
abundance distribution (full line) and with the homogeneous (log(5r/A^totai) = -8.00, log(fia/A^totai) = -9.00) abundances (dashed 
line). 



